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ABSTRACT 

The 11.2 pm emission band belongs to the family of the ‘Unidentified’ Infrared (UIR) 
emission bands seen in many astronomical environments. In this work we present a 
theoretical interpretation of the band characteristics and profile variation for a num¬ 
ber of astrophysical sources in which the carriers are subject to a range of physical 
conditions. The results of Density Functional Theory (DFT) calculations for the solo 
out-of-plane (OOP) vibrational bending modes of large polycyclic aromatic hydrocar¬ 
bon (PAH) molecules are used as input for a detailed emission model which includes 
the temperature and mass dependence of PAH band wavelength, and a PAH mass 
distribution that varies with object. Comparison of the model with astronomical spec¬ 
tra indicates that the 11.2 pm band asymmetry and profile variation can be explained 
principally in terms of the mass distribution of neutral PAHs with a small contribution 
from anharmonic effects. 
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1 INTRODUCTION 

The discovery by Gillett, Forrest & Merrill (1973) of an emis¬ 
sion feature near 11.2 pm in spectra of planetary nebulae 
opened a new era in the study of interstellar matter. It is 
now recognised as one of the strongest of the ‘unidentified 
infrared’ (UIR) or ‘aromatic infrared’ (AIB) bands that lie 
between 3 and 20 pm and which are generally attributed 
to vibrational transitions of polycyclic aromatic hydrocar¬ 
bon (PAH) molecules (for a review see Tielens 2008, and 
references therein). However, in no case has it proved possi¬ 
ble to identify a specific PAH molecule. This is unfortunate 
because knowledge of PAH shape, size distribution and de¬ 
gree of hydrogenation and/or ionisation would allow much 
greater exploitation of the spectra as probes of astrophys¬ 
ical conditions and processes. The advent of infrared (IR) 
satellites including ISO, Spitzer (Werner et al. 2004a) and 
AKARI (Murakami et al. 2007) has revealed the richness of 
the AIB spectrum and also significant differences in the pro¬ 
files of individual bands ( e.g. van Diedenhoven et al. 2004; 
Berne et al. 2007; Rosenberg et al. 2011, 2012; Boersma, Ru¬ 
bin & Allamandola 2012; Boersma, Bregman & Allamandola 
2013). Interpretation of these variations in terms of the phys¬ 
ical and chemical properties of PAHs and the astronomical 
objects in which they are found could assist in the identifica- 
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tion of PAH sub-groups, and thus narrow down the number 
of possible PAHs as carriers of AIBs. 

In this paper we focus on the 11.2 pm feature, which is 
one of the most distinctive bands in the AIB spectrum. It 
has been assigned to the C-H out-of-plane (OOP) bending 
mode of solo-containing neutral PAHs (Hony et al. 2001; 
van Diedenhoven et al. 2004; Bauschlicher, Peeters & Al¬ 
lamandola 2008; Ricca et al. 2012), although their precise 
molecular shapes and size distribution has not been estab¬ 
lished. Weaker emission features around 11.0 pm are gener¬ 
ally attributed to the same type of transition in PAH cations 
(Sloan et al. 1999; Hudgins & Allamandola 1999; Boersma, 
Bregman & Allamandola 2013). We present here an emission 
model based on Density Functional Theory (DFT) calcula¬ 
tions of vibrational transitions and intensities for a set of 
PAH molecules. The model considers the emission process 
following optical/UV excitation using the relevant stellar 
spectral energy distribution (SED) and includes the tem¬ 
perature dependence of the emission wavelength (band po¬ 
sition) as the PAH molecule cools through emission of IR 
photons. The results are used to explore the variation of 
the 11.2 pm profile in a range of astronomical objects with 
particular reference to the influence of the PAH mass distri¬ 
bution. 

The paper is arranged as follows. In Section 2 the char¬ 
acteristics of the 11.2 pm band are reviewed. The theoretical 
approach and the results of a series of DFT calculations are 
presented in Section 3 and the emission model is described in 
Section 4. Sections 5 and 6 contain the results and discussion 
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of their astrophysical implications. A possible contribution 
from acenes to the 11.0 pm feature in the Red Rectangle is 
also discussed. 

2 CHARACTERISTICS OF THE 11.2 pm BAND 

2.1 Classification of the 11.2 pm band 

The 11.2 pm feature has an asymmetric shape with a steep 
blue side and an extended red tail (Roche, Aitken & Smith 
1989; Witteborn et al. 1989). Following analysis of ISO ob¬ 
servations of objects in the Galaxy, van Diedenhoven et al. 
(2004) proposed two classes of 11.2 pm sources (An .2 and 
B 11 . 2 ) and introduced a description using a tail to top ra 
tio, i.e. the strength of the red wing relative to the peak 
intensity of the overall profile with maximum near 11.2 pm. 
Sloan et al. (2007) introduced a third class, C 11 . 2 . The class 
characteristics, illustrated with examples in Fig. 1, are: 

An .2 - the most common, with a peak wavelength ranging 
between 11.20 and 11.24 pm and a relatively short red tail. 
Sources belonging to this class have a low value of tail-to- 
top ratio and comprise various types of interstellar matter 
- reflection nebulae, H 11 regions and the general interstellar 
medium. 

B 11.2 - less common, with a peak wavelength of 11.25 pm 
and a long tail. Objects in this class, such as the Red Rect¬ 
angle (HD 44179), have a high value of tail-to-top ratio and 
are usually associated with circumstellar matter. 

C 11.2 - only a few examples, with peak wavelength rang¬ 
ing between 11.35 and 11.40 pm and an approximately sym¬ 
metric shape. Bands belonging to this class are found in 
carbon-rich objects. 

Some objects (mostly planetary nebulae) are classified as 
class A(B) 11.2 and have a mixture of the characteristics 
of the A and B classes, with a peak wavelength of class 
An .2 and a tail-to-top ratio of class Bn .2 objects (see for 
example BD+30° 3639 in Fig. 1). Also, within the single 
reflection nebula NGC 7023, variation in the 11.2 pm spectra 
from type A to A(B) has been found (Boersma, Bregman 
& Allamandola 2013). The A/B-type classification has also 
been applied to other spectral features in the 3 pm and 6 - 
9 pm regions. 

Based on a study of mid-IR Spitzer spectra of carbon- 
rich post-AGB stars in the LMC, Matsuura et al. (2014) 
have introduced separate classifications for the 6-9 pm and 
10-14 pm regions, where the spectra for the latter region are 
classified as a, (3, 7 and 5. As the spectra we consider here 
are high-resolution data from ISO-SWS (de Graauw et al. 
1996), we have elected to use the An. 2 , Bn. 2 and A(B) 11.2 
classification of van Diedenhoven et al. (2004). 

2.2 Asymmetry and peak wavelength of the 

11.2 pm band 

Interpretation of the shape and the variation in the profile of 
the 11.2 pm band is a significant challenge which requires an 
explanation of (a) the steep short-wavelength side, (b) the 
variable peak wavelength and (c) the extension to longer 
wavelength. Considering these in turn: 

(a) The wavelength at half peak intensity (A 1 / 2 ) on the 


short-wavelength side of the 11.2 pm feature in 16 spectra of 
class An. 2 sources is almost invariant, which we determine 
to be 11.171 (2a = 0.002) pm; the origin of this striking 
property is not understood. If the short-wavelength edge 
is considered to be part of a Lorentzian profile this would 
correspond to a FWHM of only c. 4 cm -1 which is narrow 
compared with other UIR bands. 

(b) The peak wavelength of the 11.2 pm feature has 
been found to depend on the effective temperature, T e ff, of 
the exciting star (Sloan et al. 2007; Keller et al. 2008). Fig. 2 
collects together data from various authors. The peak wave¬ 
length of the 11.2 pm band is independent of temperature for 
class An .2 objects (squares), but starting from class Bn .2 ( 
triangle - The Red Rectangle) it shifts to longer wavelength 
with decrease in stellar temperature. This change has been 
interpreted in terms of relatively newly formed material near 
to low temperature carbon-rich or Herbig Ae/Be stars (class 
C 11 . 2 ) compared with older processed material in the higher 
temperature environments of planetary nebulae and H 11 re¬ 
gions (Sloan et al. 2007). 

(c) Barker, Allamandola & Tielens (1987) suggested 
that the long red tail of the 11.2 pm band may be due to 
superimposed transitions between higher vibrational levels 
such as v = 2 —>- 1. In the absence of laboratory data this 
contribution was incorporated assuming that the v = 2 —> 1 
band falls at an arbitrarily chosen wavenumber 5 cm -1 
lower than for the fundamental transition. More recently, 
van Diedenhoven et al. (2004) noted that while a feature 
of class A(B)n .2 could be considered as coming from the 
same PAH population distribution as that which gives rise 
to class An .2 - but with higher internal energy, they com¬ 
mented that the difference between the An .2 and Bn .2 
classes could not readily be explained by anharmonic effects 
alone. A contribution to emission in the red tail could poten¬ 
tially come from duo-Hs of compact PAH cations (Hudgins 
& Allamandola 1999), C-H OOP transitions in PAH anions 
(Bauschlicher, Peeters & Allamandola 2009) or from pro- 
tonated PAHs (Knorke et al. 2009). In addition Wada et 
al. (2003) have discussed the possible role of 13 C isotopic 
substitution. Blind spectral decomposition studies indicate 
that a contribution to the red wing of ‘ 11.2 pm’ emission 
arises in carriers well removed from the exciting star in e.g. 
NGC 7023; this emission has been ascribed to ‘Very Small 
Grains’ (VSGs) (Rapacioli et al. 2006; Berne et al. 2007; 
Rosenberg et al. 2011, 2012). 

The shapes of the AIBs have been modelled (Schutte, 
Tielens & Allamandola 1993; Cook & Saykally 1998; Ver- 
straete et al. 2001; Pathak & Rastogi 2008); the most de¬ 
tailed study of the 11.2 pm band is that of Pech, Joblin 
& Boissel (2002) and includes both intermode anharmonic- 
ity (which from laboratory measurements of Joblin et al. 
(1995) allows the effect of cooling on emission wavelength 
and linewidth to be estimated), and intramode anharmonic- 
ity which is a contribution from vibrational transitions of 
the type v —>■ v — 1 with v ^ 2. A mass distribution for the 
PAHs was included of the form A ^) 3 ' 5 where Nc is the num¬ 
ber of carbon atoms. This affects the response of the PAHs 
to optical excitation and their subsequent cooling. However, 
in considering molecular diversity, it was assumed that all 
PAHs contributing to the 11.2 pm band profile had exactly 
the same infrared active mode frequency. It was deduced 
that the asymmetric profile appeared to be characteristic of 
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Figure 1. Comparison of the 11.2 (im profiles of class An .2 
(Orion H2S1) - in red (thin line), Bn .2 (HD 44179 - the Red 
Rectangle) and A(B)n .2 (NGC 7027). These ISO spectra are 
continuum-subtracted and normalised to the peak intensity. 
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Figure 2. Peak wavelength (Ac) and class-type of the 11.2 pm 
feature vs. effective temperature (T e //) of the central star. 
The squares are a sample of type A/A(B)n .2 sources from 
Hony et al. (2001): GGD 27-ILL (star formation region), CD- 
42 11721 (Herbig Ae/Be star), NGC 7023 (reflection nebula) and 
IRAS 21282+5050 (planetary nebula). The triangles are The Red 
Rectangle (RR) and mostly carbon stars (class Bn .2 to C 11 . 2 ) 
from Sloan et al. (2007), the 5-points stars represent Herbig 
Ae/Be stars (class Bn .2 to B/C 11 . 2 ) from Keller et al. (2008) and 
the circles are proto-planetary nebulae (class B/C 11.2 to C 11 . 2 ) 
from Zhang, Kwok & Hrivnak (2010). 


the anharmonicity of molecular modes, and that including 
a spread of frequencies according to molecular size and ge¬ 
ometry would lead to a change in overall band shape and 
consequent difficulty in achieving a good fit to the observed 
band. This contrasts with the work reported here where we 
find that the dependence of the vibrational transition fre¬ 
quencies on the masses of the 11.2 gm carriers and their rel¬ 
ative abundance are key factors in determining the overall 
11.2 fi m band profile and its variation between objects. 

3 DFT CALCULATIONS 
3.1 Theoretical Approach 

A large number of DFT calculations of PAH vibrational fre¬ 
quencies have been reported. Attention has generally been 
focused on neutral and singly ionised (i.e. radical cation) 
PAH molecules, with extension to hydrogenated, proto- 
nated, irregular and other closed-shell charged PAHs (Hud¬ 
gins, Bauschlicher & Allamandola 2001; Beegle, Wdowiak & 
Harrison 2001; Pathak & Rastogi 2005, 2006; Bauschlicher, 
Peeters & Allamandola 2008, 2009; Bauschlicher et al. 2010; 
Ricks, Douberly & Duncan 2009; Hammonds, Candian & 
Sarre 2011; Ricca, Bauschlicher & Allamandola 2011; Ricca 
et al. 2012; Boersma et al. 2014; Candian, Sarre & Tielens 
2014). The work reported here is in two closely related parts: 

First, we report the results of DFT calculations on a 
set of medium-sized PAHs, inspired in part by the work of 
Pathak & Rastogi (2005) on a series of acenes up to hep- 
tacene (see their fig. 7 which shows convergence of the solo 
OOP wavelength at high mass). We have undertaken a sys¬ 
tematic investigation to explore how the frequency and in¬ 
tensity of the out-of-plane (OOP) bending solo modes of 
neutral PAHs vary with PAH size and shape. The calcula¬ 


tions were performed using Q-Chem (Shao et al. 2006) with 
the B3LYP functional (Becke 1993; Stephens et al. 1994) 
and the 6-31G* basis set. Secondly, DFT calculations were 
carried out for selected highly symmetric compact large solo- 
containing PAHs and these results applied in the full spectral 
profile modelling. Gas-phase laboratory data were used to 
obtain mode-specific vibrational scaling factors. 

In the first part of the study PAH structures are labelled 
(x, y) where x is the number of 6 -membered rings in a row 
and y is the number of rows. Examples include anthracene 
(C 14 H 10 ) for which x = 3 and y = 1 and anthanthrene 
(C 22 H 12 ) for which x = 3 and y = 2 (see Fig. 3(a) and 
Tab. 1(a)). Vibrational frequencies computed in DFT calcu¬ 
lations are systematically high and a scaling factor is often 
invoked to bring calculation and experiment into agreement. 
Here a scaling factor for the specific solo OOP mode was em¬ 
ployed. For y = 1, 2 and 3 PAHs a scaling factor of 0.987 
was adopted in the modelling, determined by comparing the 
high-resolution gas-phase experimental value for anthracene 
of 876.7 cm -1 (Cane et al. 1997) with the unsealed DFT 
result of 888.2 cm -1 . 

The second part of our study focussed on large compact 
PAHs (see Fig. 3(b) and Tab. 1(b)). These molecules were 
chosen to investigate whether the OOP frequency shares a 
similar mass-dependent behaviour as for low-medium mass 
PAHs. The 11.2 pm band is generally attributed to large 
neutral PAHs (Hony et al. 2001; Bauschlicher, Peeters & 
Allamandola 2008; Ricca et al. 2012), their compactness 
and high symmetry making them optimal candidates to sur¬ 
vive in harsh astronomical environments. The study was 
again undertaken with B3LYP/6-31G*, using Gaussian 03 
(Frisch et al. 2004). A scaling factor of 0.975 was used for 
these molecules, deduced by reference to the experimental 
gas-phase data for the solo OOP bending mode of ovalene 
(C 32 H 14 ) (Joblin et al. 1994). 
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C66H20 



Figure 3. Molecular structures (a) for PAHs with x = 3-4 and 
y = 1-3, where x is the number of six-membered rings in a row 
and y is the number of rows. For y = 3, x indicates the number 
of edge rings in the top and bottom rows. In (b) four large PAH 
structures with D 2 h and Dg h symmetry are shown. 


The use of a scaling factor is normal practice when com¬ 
paring laboratory and DFT results. Typically, a scaling fac¬ 
tor of 0.958 for B3LYP/4-31G is used for simplicity for the 
entire spectrum, based on comparison with matrix-isolation 
data. However, it is known that this technique introduces 
unpredictable shifts in the band positions (Langhoff 1996). 
In contrast to previous studies, the scaling factor employed 
in this work is referenced to gas-phase data and is specific to 
the vibrational mode (solo OOP transitions of PAHs). For 
our analysis, this is preferable to the use of a generic scaling 
factor. At the present time there is insufficient laboratory 
gas-phase data on a range of PAHs to evaluate the appli¬ 
cability of the same scaling factor for all PAHs. However, 
the use of the well characterised B3LYP/6-31G* and the 
mode-specific scaling factor should reduce the uncertainty 
significantly. 


3.2 Results & Discussion 

Tab. 1 collects together the computed intensities and un¬ 
sealed vibrational wavenumbers of the studied molecules. 
Fig. 4 shows that for the acenes (x = 3 — 8 , y = 1) the wave¬ 
length decreases with increase in x (or number of carbon 
atoms), converging towards a high-mass asymptotic value 
of 11.03 pm. A similar characteristic holds for zig-zag edge 
molecules with two rows (stars) and three rows (filled trian¬ 
gles), but they reach the high-mass limiting value at longer 
wavelength (Fig. 4, lower panel). In all cases, the wavelength 
for the solo OOP bending mode moves to shorter wavelength 
as the number of carbon atoms (i.e. the mass) increases (see 
Tab. 1 (a) and (b)). A similar dependence of the CH OOP 


Table 1. Computed scaled DFT (B3LYP/6-31G*) C-H out-of- 
plane (OOP) vibrational wavenumbers for (a) acenes with x = 3- 
8, y = 1 and multi-row PAHs with x = 3-7, y = 2 and 3 (sf=0.987, 
see text), and for (b) large compact PAHs (sf=0.975, see text) 
shown in Fig. 3. 


(x,y) 

Molecule 

Wavenumber 

(cm -1 ) 

Wavelength 

(pm) 

Intensity 
(km mol -1 ) 

(a) 

(3,1) 

c 14 h 10 

876.7 

11.410 

55.0 

(4,1) 

C18H12 

895.3 

11.169 

69.3 

(5,1) 

C22H14 

901.7 

11.090 

87.8 

(6,1) 

C26H16 

905.7 

11.041 

106.7 

(7,1) 

C30H18 

906.6 

11.030 

124.8 

(8,1) 

C34H20 

906.6 

11.030 

146.2 

(3,2) 

C22H12 

875.2 

11.426 

96.9 

(4,2) 

C28H14 

881.8 

11.340 

109.3 

(5,2) 

C34H16 

886.5 

11.280 

123.4 

(6,2) 

C40H18 

888.1 

11.300 

133.7 

(7,2) 

C46H20 

888.1 

11.300 

172.0 

(3,3) 

C32H14 

884.2 

11.310 

123.0 

(4,3) 

C40H16 

890.5 

11.230 

155.4 

(5,3) 

C48H18 

891.3 

11.220 

185.0 

(6,3) 

C56H20 

891.3 

11.220 

204.6 

(7,3) 

C64H22 

891.3 

11.220 

231.8 

(b) 


C32H14 

877.0 

11.403 

120.3 


C42H16 

882.2 

11.328 

154.3 


C54H18 

885.0 

11.300 

190.6 


C66H20 

887.9 

11.262 

203.61 


C80H22 

889.3 

11.245 

234.31 


C96H24 

891.1 

11.222 

265.16 


frequency on PAH mass is found in the computed data of 
Bauschlicher, Peeters & Allamandola (2008) who used the 
B3LYP functional and a (smaller) 4-31G basis set; the val¬ 
ues for C 54 H 18 and C 96 H 24 are 11.061 and 10.959 pm, re¬ 
spectively. 

Larger PAH molecules with Nc > 100 are needed 
to determine the value at the asymptotic limit but this 
was beyond our computational capacity. Hence a functional 
form was used to fit the available data (once scaled) from 
Tab. 1(b) and so deduce the asymptotic wavelength limit 
for the largest PAHs: u = (b + a/Nc 2 ) -1 , where uj is the 
mode wavenumber and Nc the number of carbon atoms. 
The parameters of the best fitting are a = 0.021 ± 0.002 
and the asymptotic value b = 11.2139 ± 0.0007 pm (or 
891.75 ± 0.06 cm -1 ). Thus, the inferred asymptotic limit for 
large PAHs falls to shorter wavelength compared with low- 
medium mass PAHs, i.e. PAHs with multiple rows with y 
= 2 and 3 as shown in Fig. 3(a). The value of 11.21 pm falls 
within the range of peak wavelengths (11.20 to 11.24 pm) 
in class An. 2 objects (van Diedenhoven et al. 2004) and 
falls very close to the observed value of A 1/2 of 11.171 pm 
discussed in Section 2.2(a). 

Key results from these calculations are a) the more ex¬ 
tended (acene) and larger PAH molecules have solo out-of¬ 
plane modes at the shortest wavelengths, b) the transition 
wavelengths for acenes and for large PAHs reach asymptotic 
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be written: 


U(T) = J2 

i=l 


he LOi 

exp(hcu}i/kBT) — 1 


where ua is the wavenumber of the ith vibrational mode 
and n is the total number of modes of the molecule. The 
PAH molecule cools through its various vibrational modes, 
through a so-called radiative cascade. For the ith mode the 
emission rate 4>i is given by: 

A 1 ’ 0 

<t>i = -i- 

exp(hcu>i/kBT) — 1 

where the Einstein coefficients A 1,0 for spontaneous emission 
can be calculated from the DFT transition intensities (Cook 
& Saykally 1998). 

The fractional energy emitted in the ith mode, corre¬ 
sponding to a fall in internal energy 5U, is: 

6Ei{T)= ?* XUi 5U ( T ) 

}^i=l <Pi X Ui 


Figure 4. Upper panel: ISO-SWS spectrum of the 11.0 and 
11.2 pni band complex in Orion H2S1, an example of class An. 2 . 
Lower panel: wavelength behaviour of the OOP bending transi¬ 
tion with respect to number of carbon atoms for y = 1 (filled 
circles - acenes), y = 2 (stars) and y = 3 (filled triangles) - see 
Fig. 3(a) and Tab. 1(a). 

values at high mass, c) increasing the number of rows (larger 
PAHs) results in a faster approach to the asymptotic limit 
than for the single-row acenes, d) the higher-mass acenes 
may contribute to 11.0 pm emission. The presence of a high- 
mass asymptote for the wavelength of the largest PAHs pro¬ 
vides an attractive explanation for the steepness and wave¬ 
length invariance of the short-wavelength side (A 1 / 2 ) among 
class An .2 objects. In the following sections we explore this 
further using a detailed theoretical model. 


4 EMISSION MODEL 

In this section we describe how the DFT results for fre¬ 
quencies and intensities are incorporated into a model for 
the AIB emission occurring in astrophysical environments. 
As far as we are aware this is the first attempt to include 
mass-dependence of the PAH emission wavelength into a 
model and using, where available, reference to experimental 
results. We have chosen to employ a model based on the 
thermal approximation (as employed by e.g. Pech, Joblin & 
Boissel 2002; Pathak & Rastogi 2008). 

4.1 Theoretical Approach 

Leger, d’Hendecourt & Defourneau (1989) proposed that 
the thermal approximation can be used to describe the in¬ 
frared cooling of PAHs where this treatment assumes that a 
molecule can be considered as a heat bath with an average 
molecular energy U and temperature T. Following absorp¬ 
tion of a UV photon, a PAH molecule has an internal energy 
U(T P ), where T p is the (initial) peak temperature. 

In the harmonic approximation, the internal energy can 


The total emitted energy is obtained by integration over 
the temperature range from T p to 50 K, and weighted by 
the rate of photon absorption 

/*13.6 r>d 

Rabs = J 

where v is the frequency of the absorbed photon, a v is 
the frequency-dependent photo-absorption cross-section and 
13.6 eV represents the high-energy cut-off in the radiation 
field. For each molecule a v was taken from the French-Italian 
database 1 (Malloci, Joblin & Mulas 2007). For a molecule 
(PAH1) for which there was no entry, the cross-section was 
estimated using data available for the molecule closest in size 
(PAH2), scaled by A<? ah1 /A£ AH2 (Mulas et al. 2006). The 
excitation source is represented by a diluted Planck function 
(see Tab. 2). At high photon energy, photoionisation of 
the neutral PAH can also be significant and leads to a re¬ 
duction in the quantum yield for IR emission. To take this 
process into account, we incorporated in our modelling the 
energy-dependent photo-yield using the empirical law given 
by Le Page, Snow & Bierbaum (2001) and ionisation energies 
(IE) from the French-Italian database. When experimental 
values were not available, a theoretical vertical IE was used. 

4.2 Experimental input 

Joblin et al. 1995 studied the influence of temperature on IR 
absorption spectra of gas-phase PAHs in the laboratory and 
found a linear dependence on temperature for both band 
position and band width, arising mostly from anharmonic 
coupling between modes. The following relations for the fre¬ 
quency cj(T), and width Aoj(T), were used in the modelling: 

uj(T) = oj 0 + A uj RS + x'T (1) 

and 

A uj(T) = Atn(0) + x'T (2) 

1 http://astrochemistry.oa-cagliari.inaf.it/database. 
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Table 2. Physical conditions for objects considered in the mod¬ 
elling. T e ff represents the effective temperature of the central star 
in K, L the luminosity in units of solar luminosity (Lq), W d u the 
geometric dilution factor and D the distance in kpc. In the case 
of the HD 44179, the two components of the binary stellar system 
are taken into account. 


Object 

Orion H2S1 

NGC 7027 

HD 44179 

ClassO) 

An .2 

A(B)n .2 

Bn .2 

T ef f (K) 

4 ■ 10 4 b 

1.6 ■ 10 s d 

8250^ 6-10 4/ 

L (L 0 ) 

1.35 ■ 10 7c 

10 e 

6 10 3 I 100^ 

w dil 

2.09 • 10” 11 

b- 

1 

O 

00 

C4 

4.8 • 10“ 6 3.7- 10" 13 

D (kpc) 

0.45 6 

i e 

0.71^ 


References: (a) van Diedenhoven et al. 2004, (b) O’Dell 2001, 
(c) Verstraete et al. 2001, (d) Beintema et al. 1996, (e) Bujarrabal 
et al. 2001, (f) Men’shchikov et al. 2002. 

Table 3. Empirical coefficients for the temperature-dependent 
wavenumber and linewidth of the solo OOP mode (see equations 
(1) and (2). 



Large' 

Acenes-*- 

X ' (cm _1 /T) 

-0.0114 

-0.0065 

A 10 RS (cm” 1 ) 

-12.91 

-9.54 

X" (cm-‘/T) 

0.0157 

- 

Aoj( 0) (cm -1 ) 

0.54 

- 


Note: t Measurements on neutral gas-phase ovalene in the 550- 
820 K range and in a Ne matrix at 4 K (Joblin et al. 1994, 1995; 
Pech, Joblin & Boissel 2002). 

1 Measurements on neutral gas-phase anthracene (Semmler, Yang 
& Crawford 1991; Califano 1962; Cane et al. 1997) and in an Ar 
matrix at 10 K (Hudgins & Sandford 1998). No information is 
available on the temperature dependence of the linewidth; a value 
for anthracene in the gas phase (at 370 K) is known (Cane et al. 
1997). 

where A lors = wi( 0 ) — too can be considered to be 
an empirical red-shift (RS) between the calculated DFT 
frequency (at zero Kelvin) cu(0) and the frequency oj_l( 0), 
at zero Kelvin inferred from the experimental (laboratory) 
studies (Cook & Saykally 1998). Tab. 3 collects the coef¬ 
ficients for the solo out-of-plane mode. We consider ova¬ 
lene and anthracene to be representative of large and small 
(acene) PAH molecules, respectively. In the case of ovalene, 
x" was derived after removing the rotational contribution 
(see Fig. 8 of Pech, Joblin & Boissel 2002). For a molecule 
such as ovalene and based on experiments in absorption, 
it is expected that the emission will move to longer wave¬ 
length (lower to) as the cooling cascade occurs, and that the 
linewidth will increase. 

4.3 Single molecule calculation: influence of 
temperature? 

We consider here the emission profile resulting from the 
cooling cascade for a single PAH molecule, circumcoronene 
(C 54 H 18 , Fig.3 b), in radiation fields with T e //= 8,250 K 
and 160,000 K corresponding to HD 44179 and NGC 7027, 
respectively. The results are given in Fig. 5 and show that 



Figure 5. Emission profile computed for C 54 H 18 exposed to stel¬ 
lar radiation with T e // — 8,250 K (dotted line) and 160,000 K 
(full line). The linewidth is taken to be constant with a FWHM 
of 5 cnr 1 . 


(a) the peak wavelength is red-shifted and (b) the redward 
side extends to longer wavelength for a PAH with a high ini¬ 
tial internal temperature (i.e. irradiation as in NGC 7027). 
However, this is not in agreement with observation. The ob¬ 
served peak wavelength in NGC 7027 is lower and the red 
wing is no more extended than in HD 44179 (see Fig. 6 ). 
We also remark that a red tail originating from ‘hot-band’ 
emission from higher vibrational levels is difficult to rec¬ 
oncile with the presence of Class An .2 features (with a low 
tail-to-top ratio) in many compact and extended H 11 regions 
(van Diedenhoven et al. 2004). 

We conclude that a factor other than the radiation tem¬ 
perature is responsible for the change in 11.2 pm profile 
between objects and, given the DFT results of Section 3., 
propose that the dependence of the emission wavelength on 
PAH mass is the largest single factor in determining the 
differences in the 11.2 pm profiles between objects. This is 
now explored through modelling of astronomical emission 
profiles which includes the wavelength dependence on PAH 
mass. 


4.4 Summary of calculations 

The procedure to calculate the energy E(T) emitted in a 
band is described in Sec. 4.1. Each E(T) is then associated 
with a Lorentzian profile so that the total energy emitted in 
the solo out-of-plane bending mode is described by: 

p r \ '_ (F/2 )E m {T) _ 

SoU) ^ m ^((r/2 r + (to-to m (T)P) 

where F = Atu m (T) is the FWHM, a m is a multiplica¬ 
tive coefficient and m is summed over the molecules con¬ 
sidered. This calculation is then weighted by R a bs and a 
y 2 —minimisation routine employed to determine the best 
fit to the astronomical ISO spectra of the three sources - 
Orion H2S1, NGC 7027 and HD 44179 (see Tab. 2) - with 
the smallest possible number of contributing molecules. The 
only variable in the y 2 — minimisation is a m which is the 
number of molecules of type m included in the fitting. 
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5 ASTRONOMICAL PROFILE MODELLING 

In this section we describe the results of modelling the 
11.2 (rin band for three profile classes: (An. 2 , A(B)n .2 and 
B 11 . 2 ). It was found that only four very symmetrical solo- 
containing molecules (C 32 H 14 , C 42 Hi 6 , Cs 4 Hi 8 and C 66 H 20 , 
Tab.l b) were needed to produce acceptable fits, the sole 
variable in the \ 2 optimisation being the relative propor¬ 
tions of these PAHs ( a m ) which are taken to be representa¬ 
tive of the PAH mass distribution. However, we emphasise 
that this does not represent identification of these particular 
four PAHs; rather the results show that solo-containing neu¬ 
tral PAHs in this absolute and relative size range provide a 
good fit to the 11.2 pm feature and its profile variation with 
object. This aspect is discussed further in Section 6.1. Two 
closely related models were developed. 


5.1 Model 1: 11.2 pm modelling using 
temperature-dependent linewidths 

In the first application of the above approach ( short-dashed 
line in Fig. 6 ), the experimentally-based temperature- 
dependence for the band frequency and width (Eq. 1. and 
2. respectively - see Section 4.2) for ovalene were taken to 
hold for the larger molecules in the set {i.e. independent 
of PAH mass). However, as each molecule has a different 
initial peak temperature (determined by the stellar SED) 
and a different vibrational solo-OOP frequency accord¬ 
ing to molecular size, the cooling cascade and resultant 
emission profile necessarily differs in each case. A \ 2 
minimisation was undertaken with the only variable being 
the relative abundance of the four PAHs. This model 
provides a reasonable match with the peak wavelength and 
long-wavelength tails of the bands, but it does not account 
very well for the steep short-wavelength side seen in many 
astronomical spectra including Orion H2S1 and NGC 7027. 
This discrepancy falls in the part of the profile where 
high-mass PAHs are expected to contribute most strongly 
and lies near to the asymptotic limit discussed in Section 3.2. 

Considering the experimental observations for the 
widths of the infrared active C-H stretch transitions in 
naphthalene, pyrene and coronene 2 (Equation 2), then for a 
given temperature the width of the band is found to decrease 
with increase in PAH mass (see Tab. 2 and Fig. 4 in Joblin et 
al. 1995); for example at 600 K the linewidth values for naph¬ 
thalene, pyrene and coronene are 40, 26 and 15 cm -1 , respec¬ 
tively. The solo-OOP transition linewidth behaves similarly; 
at 370 K the experimental solo-OOP linewidths in absorp¬ 
tion for anthracene (Cane et al. 1997) and ovalene (Joblin, 
C., personal communication) are 25 cm -1 and 6 cm -1 , re¬ 
spectively. Additionally, as higher mass PAHs do not reach 
as high an internal temperature - for the same exciting SED, 
high-mass PAHs have narrower absorption (and emission) 
band widths than their lower mass counterparts. This as¬ 
pect is addressed in Model 2. 


2 Ovalene was not considered because the law was deduced from 
only two experimental points. 



10.8 11.0 11.2 11.4 11.6 11.8 

Wavelength (^m) 


Figure 6. Continuum subtracted ISO-SWS spectra of objects 
of class An .2 - Orion H2S1, A(B)n .2 - NGC 7027 and Bn .2 - 
HD 44179, the Red Rectangle (grey line), superimposed with the 
emission model results using (Model 1) temperature-dependent 
Aoj(T) (short dashed line) and (Model 2) molecule-specific con¬ 
strained widths Aoj (long dashed line) - see text for details. The 
lower panel shows the (scaled) wavelength behaviour of the solo 
OOP mode of large PAHs from Table.l b. 

5.2 Model 2: molecule-dependent 

temperature-independent linewidths 

Assuming that a Lorentzian profile describes the short wave¬ 
length side of the 11.2 pm band for Class An .2 and A(B)n .2 
objects, then from profile fitting of the short-wavelength part 
of the profile we obtain a Lorentzian peak wavelength of 
11.198 pm (or 893.02 cm -1 ) and a FWHM of only 0.05 pm 
(or 4 cm -1 ). This FWHM represents the upper limit on the 
width of a single contribution to the shorter wavelength side 
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of the 11.2 (Jin band. A similar value has been suggested by 
Cami (2011). Using this information, Fig. 6 shows the best 
X 2 ~ minimised fit with constrained FWHM values of 10, 
7.5, 5 and 4 cm -1 for C 32 H 14 , C 42 H 16 , C 54 H 18 and C 66 H 20 , 
respectively. The value of 10 cnU 1 was taken for ovalene 
with reference to the laboratory absorption data (at 600 K) 
and the 4 cm -1 value as constrained by the observational 
data (see above and Section 5.1). The fit to the steep short- 
wavelength side of the 11.2 pin band for the class An .2 and 
A(B)n ,2 objects is improved, while the long-wavelength tail 
does not change significantly. Thus an accurate description 
of the shorter wavelength side of the 11.2 pm band requires 
lower linewidths than experimentally found for the FWHM 
of e.g. ovalene. We also note the Red Rectangle, and in part 
NGC 7027, has a discrepancy in the 11.45-11.65 pm region 
in both models. We return to this in Section 6.2. 


5.3 Effect of SED on computed 11.2 pm profiles 

In section 4.3 it was shown that for a single PAH molecule, 
C 54 H 18 , the calculated emission profile with a low T e // SED 
is narrower than for higher T e ff- To explore the effect of 
the excitation SED on computed 11.2 pm profiles, we have 
taken the (fixed) mass distribution which was determined for 
NGC 7027 using the SED for NGC 7027, and calculated the 
profile using a lower T e // SED - that of the Red Rectangle. 
The calculated result in Figure 7. (upper panel, full line) is 
qualitatively similar to the single-molecule calculation and 
gives a very poor match to the astronomical spectrum for 
NGC 7027 when compared with Figure 6 . middle panel, 
Model 2. Similarly, taking the (fixed) mass distribution as 
deduced for the Red Rectangle, using a liigh-temperature 
SED (that of NGC 7027) also yields a poor match (Figure 
7. lower panel, full line). Clearly the SED used in the calcu¬ 
lation plays a role in determining a computed 11.2 pm pro¬ 
file but the SED is fixed for a given object and the observed 
profile is determined principally by the solo-containing PAH 
mass distribution. 

6 ANALYSIS AND DISCUSSION 

6.1 PAH mass distribution and the 11.2 pm profile 

Taking the sum of the contributions of the two lower mass 
PAHs (C 32 H 14 and C 42 H 16 ) to the 11.2 pm band fit in Model 
2 , and similarly taking C 54 H 18 and C 66 H 20 as representa¬ 
tive of higher mass PAHs, an indicative low massihigh mass 
abundance ratio can be determined with the break falling 
at Nc ~50. Using the x 2 fitted results, this ratio is 1.8:1 for 
Orion H2S1, 2.4:1 for NGC 7027 and 3.7:1 for the Red Rect¬ 
angle. From this we infer that the Red Rectangle PAH pop¬ 
ulation comprises lower mass PAHs than for NGC 7027 and 
Orion H2S1. We deduce that the shape and peak wavelength 
of the 11.2 pm feature can be described principally in terms 
of the PAH mass distribution, where high-mass PAHs are re¬ 
sponsible for the steep short wavelength side. The low vari¬ 
ance of A 1/2 and the steep blue side in An .2 objects arise as 
a natural consequence of the asymptotic limit reached by the 
solo OOP mode in high-mass PAHs. In class Bn .2 objects 
such as the Red Rectangle, the less steep short-wavelength 
side, the more extended redward tail and a band maximum 



10.8 11.0 11.2 11.4 11.6 11.8 


Wavelength (pm) 

Figure 7. Upper Panel. Calculated results for the 11.2 pm in 
the spectrum of NGC 7027 (grey) using the mass distribution in¬ 
ferred for this object from Model 2 and using the SED for the Red 
Rectangle. Lower Panel. Calculated results for the astronomical 

11.2 pm band in the Red Rectangle (grey) using the mass distri¬ 
bution inferred for this object from Model 2 and using the SED 
of NGC 7027. 

at longer wavelength suggests that the PAH mass distribu¬ 
tion is skewed towards lower masses than in class An .2 ob¬ 
jects. The class A(B)n .2 object NGC 7027 then falls between 
the class A and B extremes with a broad spread of PAH 
masses as suggested by the intermediate value of the mass 
distribution indicator introduced above. Using Spitzer ob¬ 
servations of NGC 7023, Boersma, Bregman & Allamandola 
(2013) fitted the PAH emission spectra as a function of off¬ 
set from the star HD 200775 using the NASA Ames PAH IR 
Spectroscopic Database (PAHdb) (Bauschlicher et al. 2010). 
The shape of the 11.2 pm feature evolves from class A 11.2 rel¬ 
atively near the star (c. 20” offset) towards Bn .2 as a denser 
region is approached at higher offset. It was concluded that 
the overall PAH spectra result from approximately equal 
proportions of large and small PAHs at Position I (a dense 
region) with larger PAHs being more prominent in the dif¬ 
fuse region (Position II). Our interpretation of the 11.2 pm 
band shape is consistent with these results, though our ap¬ 
proach differs in that our inferred mass distributions hold 
specifically for neutral solo-containing PAHs. 

6.2 Emission in the 11.4-11.7 pm region 

The Red Rectangle and to a lesser extent Orion H2S1 and 
NGC 7027, have excess emission relative to the models in the 
11.4-11.7 pm wavelength region. Based on the results of DFT 
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Figure 8. Experimental gas-phase absorption spectra 
(T~ 500 K) of anthracene (C 14 H 10 , top panel), 1 , 2 ,3,4- 

tetrahydroanthracene (C 14 H 14 , middle panel) and 1,2,3,4,5, 6 ,7, 8 - 
octahydroanthracene (C 14 H 18 , bottom panel). The data are 
taken from NIST MS Data Center / S. E. Stein (2011). As the 
number of attached hydrogens increases, the peak position of the 
solo C-H out-of-plane bending mode moves to longer wavelength. 


calculations and the gas-phase experimental OOP value for 
anthracene of 11.4 pm (at 373 K, see Fig. 6 ), it is suggested 
the excess in the Red Rectangle is due to emission from rel¬ 
atively low-mass PAHs which were not included in the mod¬ 
elling of the 11.2 pm feature presented here. The scaled DFT 
out-of-plane bending modes for anthracene C 14 H 10 (x = 3; 
y = 1 ) and anthanthrene C 22 H 12 (x = 3; y = 2 ) fall at 

11.4 pm. This too is consistent with a PAH population of 
lower mass for the Red Rectangle deduced separately in 
section 6.1. Noting the usual redward shift that occurs in 
emission, these and similar molecules probably contribute 
to emission in the 11.4-11.7 pm region. 

A second possible contribution to the ‘excess’ could 
arise from hydrogenated PAHs. Addition of hydrogen atoms 
to PAH molecules (H-PAHs) perturbs the C-H out-of-plane 
bending mode frequency, thus affecting the peak position. 
Fig. 8 shows gas-phase absorption spectra in the 11.4 pm 
region, recorded at ~500 K (NIST MS Data Center / 
S. E. Stein 2011) for the closely related molecules: an¬ 
thracene (C 14 H 10 ), 1,2,3,4-tetrahydroanthracene (C 14 H 14 ) 
and 1,2,3,4,5,6,7,8-octahydroanthracene (CmHis). Starting 
with anthracene (upper panel), increasing the number of 
hydrogen atoms (middle and lower panels) moves the peak 


position for the solo C-H out-of-plane bending mode (Hs in 
the 9,10 positions) to longer wavelength. In a major study of 
infrared spectra of hydrogenated PAHs in an argon matrix it 
was found that aromatic C-H stretching bands near 3.3 pm 
weaken and are replaced with stronger aliphatic bands near 

3.4 pm, and that aromatic C-H out-of-plane bending mode 
bands in the 11-15 pm region shift and weaken (Sandford, 
Bernstein & Materese 2013). However, very few PAHs in the 
sample set had solo-hydrogens. We tentatively suggest that 
a signal identified by Rosenberg et al. (2011) could be due 
to solo-containing hydrogenated PAH molecules rather than 
‘Very Small Grains’ - VSGs (see their Fig. 6 ); further assess¬ 
ment of this will be described elsewhere. Significantly, ob¬ 
jects classified as A(B)n .2 and Bn. 2 , show aliphatic bands at 

3.4 and 3.52 pm (Geballe et al. 1985) as would be expected 
to appear when PAH hydrogenation levels are high; these 
bands appear in the experimental spectra of hydrogenated 
anthracene, but emission in the 6-9 pm region remains weak 
as for other neutral PAHs. 


6.3 The 11-15 pm and wider spectral region 

The molecules considered here give a strong contribution to 
the UIR bands in the 11-15 pm range (Schutte, Tielens & Al- 
lamandola 1993). Fig. 9 shows the fingerprint region for the 
three prototype objects compared with the emission model 
results. Due to the lack of T-dependent experimental data 
for the position of and width of the bands in the fingerprint 
region, a single scaling factor of 0.979 and a fixed FWHM of 
12 cm -1 were applied for all bands. To combine the fit for 
the 11 . 0 - 11.2 pm band region and the remaining bands in 
the fingerprint region, we assumed that the molecular con¬ 
tribution for the solo mode in a PAH can be modelled as a 
Lorentzian function peaking at the scaled (sf= 0.979) DFT 
frequency but with a FWHM of 10 cm -1 . The predictions 
for the bands in the fingerprint region are consistent with 
the astronomical spectra of the three objects, with no more 
than a small contribution to the 12.7 pm feature. The model 
slightly overestimates the flux of the 12.0 pm band for the 
Red Rectangle, which is generally ascribed to duo hydrogen 
modes (Hony et al. 2001). However, in SWS-ISO data, the 
region around 12 pm is sensitive to the way spectra are com¬ 
bined, thus introducing extra uncertainty of 20-30 percent 
(Hony et al. 2001). A further issue is placement of the con¬ 
tinuum in the subtraction applied to the spectra. The lack of 
a detailed treatment of anharmonicity (temperature depen¬ 
dence) may also affect this modelled spectrum. It is notable 
that the model built on treatment of the 11.2 pm region does 
not introduce features which are not present in astrophysi¬ 
cal spectra of the 11-15 pm region; this is also the case for 
the 6-11 pm region. For example, the strongest calculated 
vibration of C 54 H 18 in the 6-9 pm region has a strength of 
50 km mob 1 (Candian 2012) compared with 191 km mob 1 
for the solo OOP mode, i.e. an intensity ratio of 0.3. The 
effect of the radiative cascade during the PAH emission will 
reduce this ratio even more. Also, the 6-9 pm region is gen¬ 
erally attributed to ionised PAHs, which show very strong 
C-C and C-H in plane stretching modes, compared to their 
neutral counterparts (Hudgins & Allamandola 1999). As for 
the 3.3 pm band, it is known that DFT calculations overes¬ 
timate the intensity of the C-H stretching modes (Langhoff 
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Figure 9. Emission model in the fingerprint region for the as¬ 
tronomical objects as in Fig. 6. A scaling factor of 0.979 and a 
fixed FWHM of 12 cm -1 were applied. To combine the fit for 
the 11.0-11.2 urn band complex and the remaining bands in the 
fingerprint region, the molecular contribution for the solo mode 
in a PAH was assumed to be a Lorentzian function peaking at the 
scaled (sf= 0.979) DFT frequency and with FWHM of 10 cm -1 . 

1996), making it difficult to directly compare our model with 
the astronomical spectra. 

6.4 Emission in the 11.0 pm region 

6-4-1 Emission in the 11-Op.m region - a summary 

Given the proximity of 11.0 pm emission to the 11.2 pm 
band, we review here whether there is significant evidence 
for a related origin. Roche, Aitken & Smith (1991) first iden¬ 
tified a faint feature on the short-wavelength side of the 
11.2 pm band which appears quite prominently in spectra of 
a number of objects (Hony et al. 2001). The most commonly 
occurring band has an approximately symmetric shape cen¬ 
tered at 11.05 pm, a rather narrow FWHM of 0.1 pm and 
an intensity which ranges between 2 and 10% of the 11.2 pm 


feature. In a few cases a somewhat higher 11.0/11.2 ratio is 
found: two Herbig Ae/Be stars, MWC 1080 (Sakon et al. 
2007) and HD 37411(Boersma et al. 2008), in the Hu re¬ 
gion IRAS 18434-0242 (Roelfsema et al. 1996) and in two 
reflection nebulae, NGC 1333 SVS3 (Sloan et al. 1999) and 
NGC 7023 (Werner et al. 2004b; Boersma, Bregman & Al- 
lamandola 2013). 

Emission features in the 11.0 pm region have been de¬ 
scribed as ‘blue outliers’ of the 11.2 pm band by Sloan et 
al. (1999). On increasing distance from the hot star SVS3 
it was found that a shorter wavelength part at ~ 10.8 pm 
disappeared first, whereas ~ 11.0 pm emission was more 
persistent. Using the results of theoretical calculations by 
Langhoff (1996), Sloan et al. (1999) considered whether the 
‘11.0 pm’ emission might be due to neutral acenes, such 
as anthracene or pentacene, but following comparison with 
IR matrix data they concluded that this was unlikely and 
PAH cations were preferred as carriers. Hudgins & Alla- 
mandola (1999) discussed this in some detail and extended 
consideration across the 11-14 pm region. The idea of a 
cationic origin gained support from DFT calculations on 
large solo-containing PAHs (Bauschlicher, Peeters & Alla- 
mandola 2008; Ricca et al. 2012) where the (scaled) solo 
wavelengths for the cations were found to fall to wave¬ 
lengths shorter than that for neutrals. Werner et al. (2004b) 
and Sakon et al. (2007) found that the c. 11.0 pm emis¬ 
sion reduces in intensity moving away from the exciting 
star in NGC 7023 and MWC 1080, respectively, a result 
for NGC 7023 recently reinforced by Spitzer observations 
(Boersma, Bregman & Allamandola 2013). 

The most commonly discussed carrier of the c. 11.0 pm 
emission carrier is PAH cations. However, Povich et al. 
(2007) observed a lack of variation in the shape and/or in¬ 
tensity of the 11.0 pm band in the photodissociation region 
M17-SW, even where the level of the ionisation in the gas 
varied with position. Moreover, in recent work on the evo¬ 
lution of the 11.2 pm and 11.0 pm features SE of the Orion 
Bright Bar (Boersma, Rubin & Allamandola 2012), the in¬ 
tegrated strength of the main 11.2 pm feature dropped by 
a factor of 5.7 between 2’.6 and 5’.7 whereas the 11.0 pm 
feature dropped by a factor of 2.7. This result was described 
as ‘counterintuitive’ in terms of the 11.0 pm emission being 
due to PAH cations but rationalised as being influenced by 
the degree of hydrogenation. In a recent Spitzer study of the 
mid-infrared structure of the massive star-formation region 
W49A further puzzles have emerged; the 11.0/11.2 ratio is 
found to be constant across the whole region and with a very 
low value of ~0.01 (Stock et al. 2014). The ratio is more 
commonly in the range 0.02-0.2 as discussed by Stock et al. 
(2014). In a recent study fitting the spectrum of NGC 7023 
with the NASA AMES PAH database, Boersma, Bregman & 
Allamandola (2013) noted that the 11.0 pm band could be 
reproduced by ionised nitrogen-containing PAHs. In sum¬ 
mary there are many questions as to the origin (s) of fea¬ 
ture^) in the 11.0 pm region. 

6-4-2 Red Rectangle emission in the 11.0\im region 

It was inferred earlier that the Red Rectangle has a popula¬ 
tion of lower-mass PAHs than for Orion H2S1 and NGC 7027 
(Section 6.1), so it is of interest to explore whether low-mass 
PAHs including acenes might contribute to the 11.0 pm re- 
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Figure 10. Comparison between modelled and astronomical 
spectrum of the Red Rectangle (grey line) for the 11.0-11.2 pm 
band complex. The fit takes into consideration the four large 
PAHs plus pentacene, hexacene and heptacene. The width for 
the contribution is fixed at 8 cm'" 1 for acenes and as in Model 2 
for large PAHs. The value of the low mass: high mass ratio of 3.7: 
1 deduced in section 6.1 is changed by less than 10% on including 
acenes in the mode. The insert shows a comparison between the 
contribution of quarto OOP bending modes in acenes and the Red 
Rectangle spectrum (grey line) aroun 13.5 pm. The contribution 
is fitted using the constraints described in Sec 6.3. 

gion in this object. It is notable that DFT calculations pre¬ 
dict an asymptotic wavelength for the longer acenes which 
falls at a shorter wavelength than the 11.2 pm band (Sec¬ 
tion 3.2). We have conducted a Model 2-type calculation 
including acenes which yields the result shown in Fig. 10. 
The fit to the 11.0 pm region is improved and the 11.2 pm 
profile is not affected, suggesting that, at least for this ob¬ 
ject, the acene series may make a significant contribution to 
11.0 pm emission. The lack of a very good fit between the 
11.0 pm feature and larger acenes may be due to the under¬ 
estimation of the frequency of the OOP solo mode transi¬ 
tion for these molecules. We are exploring this further by 
undertaking DFT calculations which include anharmonicity 
explicitly. 

The IR. vibrational spectrum of neutral, long acenes 
is dominated by strong C-H oop modes (solo and quarto) 
(Pathak & Rastogi 2005). We therefore checked the contri¬ 
bution of the acene quarto modes by modelling them (Model 
2) and compared the result with the emission in the Red 
Rectangle (Fig. 10, small insert). Clearly, the prediction is 
consistent with the astronomical spectrum. The overall re¬ 
sult that the Red Rectangle has a population of relatively 
low mass PAHs may be of significance in determining the 
origin of the strong unidentified optical emission bands seen 
in this object (Warren-Smith et al. 1981). 


7 CONCLUSIONS 

Based on the results of DFT calculations, we created a de¬ 
tailed emission model and applied it to a small set of PAH 
molecules with the aim of understanding profile and profile 
variations in the 11.2 pm band. 

The model shows that the 11.2 pm feature, regardless 
of the A/B classification, can be fitted by a mass distribu¬ 
tion of a few, large neutral PAHs containing solo hydrogens. 


We infer that the mass distribution is principally respon¬ 
sible for the red tail of the feature and the steep short- 
wavelength side, and hence offers a good explanation for 
the 11.2 pm band asymmetry and its variation between and 
within objects. In this picture, large neutral solo-containing 
PAHs dominate class An .2 objects and smaller PAHs class 
Bn .2 objects, whereas class A(B)n .2 possesses a more even 
mass distribution. We also note that a PAH mass range of 
30-70 carbon atoms compares favourably with the masses of 
fullerenes, recently discovered in several objects ( e.g. Cami 
et al. 2010; Sellgren et al. 2010). The present model gives 
a more detailed description of the different physical effects 
(anharmonicity, photo-ionisation and astronomical environ¬ 
ment) involved in building up the feature. 

In this paper we described the results of a detailed emis¬ 
sion model applied to the solo mode in PAHs and compared 
them to astronomical data. Based on the results we find 
that the profile of the 11.2 pm band and its variation are 
well described principally in terms of the mass distribution 
of neutral solo-containing PAHs. 
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